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Problem  Statement 

This  document  describes  ongoing  research  for  the  formulation,  analysis  and 
implementation  of  a  programmable  device,  termed  Quantum  Function  Evaluator  (QFE), 
that  uses  quantum  state  propagation  as  a  paradigm  for  computing  objects  of  a  class  I  of 
functions  herein  referred  to  as  the  computable  class.  The  proposed  research  is  divided 
into  three  phases:  Formulation,  Evaluation  and  Implementation.  We  will  discuss  an 
outline  of  the  first  phase  in  the  next  sections. 

Outline  of  formulation  Phase 

Over  a  period  of  several  years,  Kohn  and  Nerode  have  conducted  joint  and  independent 
research  to  explore  formd  methods  in  Hybrid  System  Theories  [1],  [2],  [3],  [4],  [5],  [6], 
[7]  to  determine  effective  computational  procedures  for  generating  control  laws  in  a 
variety  of  application  domains.  The  quantum  computing  paradigm  and  its  possible 
implementation  proposed  for  this  research  effort,  constitute  a  significant  extension  of  the 
results  in  our  early  efforts.  This  outline  provides  a  brief  synopsis  of  our  proposed 
paradigm  and  a  preliminary  sketch  of  a  possible  implementation  architecture.  The 
Formulation  phase  of  our  study  is  composed  of  two  major  tasks:  the  detailed  specification 
of  the  proposed  paradigm  and  the  detailed  specification  of  an  implementation 
architecture .  We  will  outline  the  major  features  of  these  two  items  next. 

Paradigm 

The  main  aspects  of  our  proposed  paradigm  are  summarized  in  Figure  1 .  The  proposed 
computational  paradigm  is  composed  of  seven  sequential  steps.  We  outline  their 
functionality  next. 

Input  Function:  Each  Input  Function  is  either  a  map  of  the  form: 

f  :  S,  X  •••  xSn  D 

where  Sj ,  i  =1, ...,  N  and  D  are  finite  subsets  of  the  naturals  of  the  form: 

S,  =  {0,  N,} 

D  =  {0,  Nj 

or  a  solution  of  an  iterative  process  of  the  form: 
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with 


yLi  =  -  y')  >  i  =  ^ 

f  :  S"  -»  S  ,  S  =  {0,1,  ...,  Ns} 

Let  p  be  the  natural  number  defined  by: 

p  =  max{{Ni ,  i  =  1,  N},  nJ 

Then,  we  can  express  f  by  the  encoding: 

F  :  [O,  1,  p''-']  ^  [0,  1,  ...,p] 

For  X  = 

s=0 

(1) 

p(x)  =  •••’  "n-i)  ’  If  defined 

0  Otherwise 


Figure  1.  Quantum  Computing  Paradigm 
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Thus,  the  encoding  is  a  discrete  function  with  domain  on  certain  points  of  the  real  line, 
taking  values  in  a  subset  of  points  of  the  real  line  coinciding  with  the  values  in  the  range 
of  f  and  mapping  to  0  any  encoding  x  that  corresponds  to  a  tuple  not  in  the  domain  of  f. 

The  central  objective  behind  our  paradigm  is  to  find  effective  and  fast  means  to  compute 
discrete  finite-valued  functions  via  Quantum  Approximations  to  their  encoding.  We  will 
characterize  the  nature  of  these  approximations  later  on.  We  conclude  the  overview  of  the 
Input  Function  element  of  our  paradigm  with  a  formal  definition  of  the  class  I  of 
computable  functions. 

The  class  I  of  computable  input  functions  is  defined  as  follows: 

(i)  I  contains  all  the  functions  encodable  in  one  step  as  above. 

(ii)  If  {f, ,  -  ,  4  ,  fi  :  S,  X  •••  xSn  ^  ,  K  finite,  S,  =  {l,  ...,  Nj,  N-,  ,  finite  }} 

is  a  subset  of  I,  so  is  their  direct  sum;  f,  ©  •  •  •  ©  f^  :  S,  x  •  •  •  x  . 

(iii)  I  contains  the  projection  functions: 

Pi  :S,  X  •••  xSn  ->Si ,  Pi  (ni ,  ...,  ni,  ...,  n^)  =ni. 

(iv) If  {f,  ,  •••  ,  4  ,  fi  :  Sj  X  •••  xSfi  S^+i,  K  finite,  S-  ={l,  ...,  Ni,  Ni  ,  finite  }} 

is  a  subset  of  I  and  g  :  is  in  I,  so  is  the  composition: 

g(f, ,  •••  ,4):  S,x  •••  xSn^S„^,. 

(v)  If  for  each  ne  N,  g :  NxS,  x  •••  xS^ -» N  is  in  I  so  is 
min{n,  g(n  ,  n, ,  •••  ,  n^)  =  O}. 

(vi) If  g:  NXSk  is  in  I,  then  the  family  of  functions, 

{f:  NxSjX  •••  xSk^  Sk+1,  f(n  +  l,  nj,...,  nj  =  g(n,f(n,n,,  ...,  nj)}  is  in  I. 

(vii)  Any  function  constructed  by  the  finite  application  of  (i)-(vi)  is  in  I. 

We  note  that  the  defining  characteristic  of  I  is  that  its  elements  are  either  an  encoding  or 
can  be  transformed  into  an  encoding  via  finite  number  of  steps  (ii)-  (vii)  above.  We 
conclude  the  outline  of  the  input  function  with  the  following  observation:  although  the 
class  I  of  computable  functions  I  includes  only  completely  specified  functions,  it  can  be 
extended  to  a  larger  class  which  include  discrete  partially  specified  functions  (relations). 
We  will  discuss  this  extension  in  a  future  report. 


Continualization:  We  can  complete  the  encoding  function  F  into  a  step  function  O  with 
the  following  identity: 


O  :[0,  p^-')  ^[O,  p'^-') 

0(y)=F(x),  for  each  y,  y  e[x,  x  +  1) ,  xe{0,  ...,p'^'} 


3 


Figure  2  illustrates  this  identity  for  the  encoding  of  the  function  given  in  Table  1.  We 

refer  to  the  process  of  constructing  O  as  Continualization .  With  some  needed 
modifications,  this  continualization  process  can  be  extended  to  discrete  functions  that  are 
specified  as  solutions  of  iterations  of  the  form; 


Table  1. 


Ho 

"1 

f 

0 

0 

2 

0 

1 

1 

0 

2 

1 

1 

0 

0 

1 

1 

2 

1 

2 

2 

2 

0 

1 

2 

1 

0 

2 

2 

0 

X 

F 

0 

2 

1 

1 

2 

1 

3 

0 

4 

2 

5 

2 

6 

1 

7 

0 

8 

0 

Original  Function 


Encoding 


Step  Function 

Figure  2.  Continualization  of  a  function 


^k+l  ~  8(^k’ 

with  (3) 

Z^.  E  S^,  S  =  {0,  N,  N  finite}  for  each  k,  k  an  integer 
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Following  Kushner  and  Clark  [8],  the  encoding  and  continualization  of  a  process  of  the 
form  of  (3)  leads  to  a  differential  equation  of  the  form  of  (4)  together  with  a  sampling 
rule  of  the  form  of  (5): 


x(t)=  G(x(t),  t)  +  B(t)  (4) 


G(-,  t)  =  continualized  version  of  g(-,  k),  t  €  [kA,  (k  +  1)A) ,  k  e  N  (5) 

and  B(t)  is  a  'convergence'  function  satisfying  the  following  conditions: 

B  is  bounded,  and  B(t)  goes  to  zero  as  t  — >  ooquadratically.  Under  general  conditions,  one 
can  show  [9]  that  particular  solutions  of  (4)  are  asymptotically  stable  in  the  sense  of 

Lyapunov.  In  (5),  A  is  a  scaling  parameter  chosen  to  transform  the  iteration  variable  k 
into  the  time  variable  t. 


Finite  Hilbert  Expansion  :  The  continualized  functions,  appearing  in  (1)  or  (4)  may  be 
expanded  in  a  discrete  Hilbert  space  in  terms  of  an  orthogonal  basis  set  such  as  the  set 

of  p-valued  Chrestenton  functions.  The  general  form  of  a  Chrestenton  basis  function  is 
given  below. 


Chrestenton  function :  ^^(x)  =  exp 


VP 


J 


For  the  sake  of  simplicity,  we  will  continue  our  discussion  here  under  the  assumption  that 
the  selected  bases  are  the  Chrestenton  functions,  although  during  the  study  of  tunneling 
as  a  mechanism  for  implementing  quantum  oracles  we  may  need  to  switch  to  other 
bases. 


Let  Bp  =  |(p’' I  (p^  :  [0,p^'']  [0,p'''‘‘] ,  k  =  0,  •••,  p^ -ij  be  an  orthogonal  basis  for 

/p.  That  is,  J^''(x)-^^(x)dx  =  c-5|j^,  for  each  k,  1  in  |0,  •••,  p^-l}.  The  function 
F  can  be  expanded  in  terms  of  Bp  as  follows: 

«>W  =  lA'p-Cx)  (6) 

xefO.p]" 

with 
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We  note  that  if  we  know  the  spectrum  [Py.]  We  can  easily  compute  O,  via  (1),  extract 
the  encoding  F  from  it  and  from  F  determine  the  original  function  f.  However  even  for 
functions  of  more  than  academic  interest,  the  number  of  terms  in  the  spectrum  runs  in  the 
range  of  10^  to  10^^  jhe  central  idea,  in  this  regard,  is  that  we  can  find  very  close 
approximations  to  this  evaluation  by  a  hardware  quantum  device,  which  we  will  discuss 
in  the  next  section.  For  the  moment,  in  this  section,  we  continue  with  the  description  of 
the  paradigm. 

Variational  Model:  Given  a  computational  process  to  be  carried  out,  expressed  in  the 
form  of  the  differential  equation  (4),  we  want  to  formulate  the  computation  of  the 
solution  x(t),  as  a  variational  problem.  Our  purpose  is  to  proceed  from  this  to  obtain  a 
quantum  mechanical  "program  "  for  computing  the  solution.  For  the  purposes  of  this 
report,  we  will  describe  the  variational  formulation  formally.  Mathematical  rigor  will 
follow  in  future  reports. 

The  variational  formulation  is  of  the  form: 


minjL(x,  x,  t)-dt  (7) 

where  L  known,  as  the  lagrangian  function,  is  chosen  so  that  for  specific  boundary 
conditions  the  solutions  of  (7)  and  (4)  coincide.  For  the  purposes  of  our  procedure  we 
assume  that  L  is  three  times  continuously  differentiable  in  its  arguments.  We  now 
proceed  to  describe  how  L  is  determined  from  the  function  G,  constructed  earlier. 
Towards  this  objective  we  write  the  necessary  conditions  for  optimality  in  (7),  The  Euler 
Lagrange  conditions.  They  are  expressed  by  the  following  second  order  differential 
equation: 


^xx  +  X  L,,  +  L„  -  L,  =  0  (8) 

Here  and  in  the  derivations  that  follow  subindices  indicate  partial  derivatives. 


Differentiating  (4)  with  respect  to  t,  to  obtain  a  definition  for  x,  replacing  in  (8),  and 
differentiating  the  result  with  respect  to  x,  we  obtain,  after  some  algebra,  the  following 
partial  differential  equation  : 

Lxxi  +  +  xL^^^  +  (xG,^  +  G,  +B,)L^^^  =  0  (9) 


Let 


q(x,  X,  t)  =  Li,(x,  X,  t) 


In  terms  of  q,  (9)  takes  the  form. 


q,  +  G,q  +  xq,  +  (xG,  +  G,  +  B,)q,  =  0 


(10) 


(11) 


That  is,  q  satisfies  a  linear  hyperbolic  equation,  whose  solution  can  be  computed  by  the 
method  of  characteristics  ([10]).  Once  a  solution  for  q  is  obtained  L  can  be  determined 
from  (10)  by  a  double  quadrature.  In  summary,  given  (4)  we  construct  by  the  procedure 
sketched  above  the  corresponding  lagrangian. 
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Quantization:  In  order  to  quantize  the  program  described  by  (7),  one  introduces  the 
^nomcal  conjugate  variables:  generalized  position  and  momentum  [12].  In  the  canonical 
Quantization  these  are  chosen  as  follows:  position  variable,  representing  the  position 
operator,  is  represented  by  x,  and  the  momentum  variable  is  chosen  as: 

P  =  L,  (12) 

The  components  of  x,  p  satisfy  the  canonical  commutation  relations: 

K,P,]  =  i^k,.  k,  1  =  1,....  ,N  (13) 

Next,  we  define  the  Hamiltonian  of  the  system  as  follows: 

H(x,  p)=  XPk-  Xk- L  (14) 

The  equation  of  motion  of  the  system  representing  the  program  to  be  executed  is  then  the 
Schodinger  equation  given  by  : 


ih^lTCt))  =  h(^x,  1|-]  !■!<(, ))  (15) 

Realization:  Assume  that  a  physical  system,  the  hardware,  is  available  to  us.  This  system 
is  characterized  by  the  Hamiltonian  operator  Hq.  The  realization  step  consists  in 

engineering  a  field  Hamiltonian  H^  such  that  the  hardware,  when  interacting  with  this 
field  satisfies  the  following  condition: 


|Ho  +  Hf  -  H|  <  e  (16) 

where  |  •  |  denotes  a  suitably  selected  operator  norm  and  e  is  an  engineering  parameter 
chosen  to  satisfy  precision  requirements  for  the  computation. 

The  general  idea  is  that  the  discrete  spectrum  of  the  composite  Hamiltonian  Hg  +  H^ 
approximates  the  spectrum  of  the  function  or  process  being  computed  (see  (6)  above). 
Precisely,  we  will  devote  our  next  report  to  prove  the  following  result: 

"Let 

k  =  1,  ...,  N}  be  a  subset  of  the  eigenvalues  of  Hg  +  Hf  with  corresponding 
eigenvectors  {^k  k  =  1,  ...,  N} : 


(  H.  +  H,  )  |¥,>  = 

then 


K  -  A|  s  ^ 

and 


1 


J|/(x)-'Pjx)|-dx  <  e 

0 

where  cp'^  is  a  Chrestenton  function" 

Thus  computing  a  discrete  function  or  a  discrete  process  can  be  approximated  by  exciting 
the  hardware  system  appropriately  and  reading  the  resulting  spectrum. 

Computation:  The  computation,  to  a  large  extent,  consists  in  simulating  the  system 
characterized  by  Hg  +  Hf  using  the  hardware  characterized  by  Hg,  and  excited  by  the 

field  whose  Hamiltonian  is  Hj.  Our  research  about  computation  is  geared  towards 
developing  an  effective  translator  that  converts  an  input  function  or  process  of  the  class  I 
into  a  Hamiltonian  operator  H  and  a  resolution  element  that  extracts  the  excitation 
Hamiltonian  Hf.  Our  long  term  objectives  are  to  design  and  implement  a  prototype 
system  that  operates  according  to  the  paradigm  outlined  in  this  section. 


Figure  3.  Functional  Architecture  for  Quantum  Computing 
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Architecture 


In  this  section  we  propose  a  preliminary  design  of  an  architecture  that  implements  the 
paradigm  presented  in  the  previous  section.  A  functional  model  of  our  proposed 
architecture  is  depicted  in  Figure  3.  We  will  next  discuss  briefly  the  functionality  of  the 
elements  composing  our  proposed  architecture: 

Input  Function  -  Given  a  discrete  function  or  process  to  be  computed,  this  element 
implements  the  compilation  process;  that  is,  it  determines  the  excitation  Hamiltonian, 
Hf,  which  encodes  the  function  or  process  .  This  is  a  symbolic  computation  whose  steps 
are  carried  out  following  the  paradigm  described  in  the  previous  section.  These  steps  can 
be  described  by  an  iterative  process  such  as  (3)  so  in  principle  the  compilation  process 
could  be  carried  out  by  the  architecture  itself.  We  will  explore  the  mechanization  of  the 
compilation  process  in  this  direction. 

Excitation  Generator  -  This  is  a  device  that  realizes  physically  the  field  excitation 
encoded  in  Hf.  That  is,  the  excitation  generator  converts  a  description  of  the 
computation  into  a  physical  implementation  .  The  idea  is  for  this  device  to  address  each 
individual  particle  of  the  quantum  process  described  below  and  to  excite  them  according 
to  the  desired  state  behavior  dictated  by: 

=  (Ho  +  H,)(x,  |T(t)> 

Our  research  in  this  area  will  be  focused  on  the  development  of  a  mechanization  process 
for  realizing  this  functionality.  We  will  consider  two  alternatives:  optical  or  particle 
(electron)  excitation.  A  more  detailed  discussion  of  these  two  alternatives  will  be 
provided  in  our  next  report. 

Quantum  Processor  -  This  is  the  device  which  actually  carries  the  computation.  Without 
excitation,  it  is  a  realization  of  the  Hamiltonian  Hq  .  After  excitation  it  is  a  realization  of 
H  =  Hq  +  Hf .  We  will  provide  next  an  abstracted  model  of  the  physical  characteristics 
of  the  Quantum  Processor. 

The  Quantum  Processor  is  an  array  of  identical  particles  assembled  into  a  regular  lattice. 
For  the  purposes  of  discussion  we  will  consider  a  two-dimensional  lattice.  A  later  report 
will  be  devoted  to  the  formulation  of  the  physical  characteristics  of  the  lattice;  in  this 
report  we  will  be  concerned  only  with  the  formulation  of  some  of  its  computational 
behavior. 

A  diagram  of  the  Quantum  Processor  is  shown  in  Figure  4.  The  device  is  composed  of 
two  elements:  the  computational  lattice  of  particles  and  the  field  excitation  device.  From 
a  computational  point  of  view  each  particle  which  is  allocated  to  a  node  in  the  lattice,  can 
be  represented  as  a  non-deterministic,  two-level  state  automaton,  and  an  interface 
function  called  Input  selector  function  as  shown  in  Figure  5.  We  proceed  to  describe  their 
functionality  next. 

In  the  Lower  Automaton,  the  block  labeled  'State  transition'  in  figure  5  characterizes  the 
programmable  discrete  spectrum  of  the  particle.  The  states  of  the  automaton  represent 
energy  levels,  and  edges  represent  allowable  energy  transitions.  This  is  illustrated  in 
Figure  6  below. 
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Boundary 

Quantum  Computational 
Element 

Quantum  forward  bond 

Quantum  lateral  bond 
Excitation  path 


Figure  4.  Quantum  Processor 


Figure  5.  Quantum  Particle 


rri 


In  Figure  6,  the  energy  levels  correspond  to  eigenvalues  ( or  sets  of  eigenvalues  forming  a 
band  grouped  together  as  a  single  eigenstate)  of  the  corresponding  Hamiltonian.  The 
edges  correspond  to  energy  transitions.  Edges  pointing  up  are  driven  by  excitation;  that 
is,  the  transition  is  effected  by  absorption  of  energy  from  the  excitation.  Edges  pointing 
down  correspond  to  relaxation  effects:  the  particle  releases  energy  of  the  appropriate 
frequency  either  to  neighborhood  particles  or  to  the  field  (see  Figure  7). 


f  1 


Energy  Level 

►  Allowable  Transition 

Figure  6.  State  energy  transition  example 


Figure  7.  A  state  transition 
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In  Figure  7,  and  ^  are  energy  levels,  and  the  transition  between  them  is  driven  by 
excitation  energy  to  frequency  Vj  2  ■  Note  from  Figure  5  that  energy  can  come  to  the 
particle  either  by  interaction  with  the  other  particles  or  from  the  external  excitation. 
Relaxation  transitions  are  similar. 

The  state  transition  in  the  lower  level  automaton  is  controlled  by  the  mixing  state 
probability  density  transition  computed  by  the  Upper  Level  Automaton.  State  transition  in 
this  automaton  is  called  'probability  transition'  in  Figure  5. 

The  mechanism  that  implements  commands  issued  by  the  upper  level  automaton  in  the 
lower  level  automaton  is  called  Transition  Selector  in  Figure  5.  We  will  explain  its 
functionality  next.  For  this  we  need  some  preliminary  definitions. 

Let  S=  ||T'^(t)),  k  e  Np,  Np  finite}  be  the  states  of  the  lower  automaton  in  each 

particle,  the  let  A  be  the  update  time  of  the  particle.  The  update  time  of  the  particle  is 
the  determined  to  be  larger  than  10  times  the  maximum  relaxation  time  of  any  of  the 
state  transitions  in  the  lower  automaton.  The  state  of  the  lower  level  automaton  in  each 

interval  [t,  t  +  A)  is  a  Chattering  Combination  of  the  elements  of  S.  A  Chattering 

Combination  of  the  set  of  functions  in  the  interval  is  a  function  |^(t))  defined  as 
follows; 


iT'(T))  = 


TsIjJt) 


where  Ij^(t),  e  Np,  is  a  semi-open  interval  in  [t,  t  + A) : 


ZA.(t),  t  +  ^A.  (t) 


(17) 


(18) 


We  note  from  (17)  that  the  composite  state  |'F(t))  is  constructed  by  'stitching'  together 

the  pure  states  in  S.  The  time  spent  in  pure  state  e  S,  AjJt)  is  a  function  of  the 
characteristic  excitation  or  relaxation  times  associated  with  the  state  [11].  We  note  that 


=  A 

(19) 

and,  if  we  define 

k 

«j(t)= 

JkV  f  ^ 

(20) 

then  equation  (19)  can  be  written  as 
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(21) 


Zai.(t)  =  1 

k 

The  set  of  all  Chattering  Combinations  of  S  in  the  interval  [t,  t  +  A)  is  denoted  by 

S(t,t  +  A).  The  union  of  these  sets  for  all  time  is  denoted  by  S.  Now  we  can  specify  the 
operation  of  the  Transition  Selector  in  the  lower  automaton  of  our  particle  model.  For 
each  interval  [t,  t  +  A),  the  Transition  Selector  receives  a  command  from  the  upper 

level  automaton,  which  consists  of  an  ordered  tuple  of  coefficients  (aj_^(t),  €Np  ^ 

satisfying  (21),  and  then  it  computes  the  mixed  state  of  the  particle,  |^(t)),  according  to 
(17).  The  state  of  the  lattice  is  composed  of  the  states  of  each  of  its  particles.  We  will  see 
shortly  that  the  Chattering  Coefficients  have  the  interpretation  of  state  occupancy 
probabilities.  To  demonstrate  that,  we  need  to  discuss  the  dynamic  structure  of  the  upper 
level  automaton. 

In  general,  there  is  not  enough  information  to  say  that  the  lattice  or  any  of  its  particles  is 
characterized  by  a  specific  state  function.  The  best  we  can  do,  in  order  to  describe  the 
computation,  is  to  give  a  probabilistic  description.  In  the  quantum  formalism,  this 
description  is  referred  to  as  the  probability  density  description  [12]. 

Let  S  be  the  set  of  primitive  states  of  the  lower  level  automaton.  Let  pj  be  the 

probability  that  the  particle  is  in  state  'Fj  e  S.  The  probability  density  operator  [13],  p, 
is  defined  by 


p(t)  =  ZPj  |'*'jW)(^j(t)|  (22) 

j 


By  differentiating  (22)  and  using  (15),  after  some  algebra,  we  obtain  an  (operator) 
differential  equation  for  p  : 


ih— p  =  Hp-pH  (23) 

The  operator  is  termed  the  commutator  and  is  written  as  [H,p].  Thus 

=  [H,  p]  (24) 

Given  an  observable  characterized  by,  say,  operator  C,  the  associated  observed  quantity, 
denoted  by  (C),  is  given  by  the  expectation  of  C  relative  to  p  : 

(C)=  trace(p-C)  (25) 

Equation  (24)  characterizes  the  computation  carried  out  by  the  upper  level  automaton. 
Equation  (25)  characterizes  the  Transducer  of  our  architecture  (see  Figure  3).  We  will 
devote  a  future  report  to  discussion  of  this  device  in  detail. 
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Notice  that  the  state  function,  computed  by  the  lower  level  automaton,  could  be  given  as 
a  linear  combination  of  the  states  in  S.  We  chose  to  model  it  as  a  Chattering 
Combination,  which  turns  out  to  be  equivalent  in  a  specific  sense  (as  we  will  show 
shortly)  because  in  this  form  it  will  allow  us  to  formulate  the  sequence  of  excitation 
steps  (realized  by  the  excitation  element,  see  Figures  4  and  5).  To  a  large  extent, 
programming  the  quantum  processor  is  tantamount  to  determining  this  sequence.  To 
justify  this  statement,  we  must  explain  the  sense  in  which  the  Chattering  and  Linear  State 
Combinations  are  equivalent,  because  an  extension  of  this  result  to  Excitation 
Hamiltonians  will  provide  us  with  a  strategy  for  implementing  excitation  sequencing  . 

The  equivalence  between  Linear  and  Chattering  Combinations  of  state  functions  from  a 
given  set  S  is  established  in  the  following  version  of  the  Chattering  Lemma  [14]: 


Chattering  Lemma.  Let  S=  ||'Fk(t)),  k  €  Np,  Np  finite}  and  let  S  be  the  set  of 
chattering  combinations  of  S.  Let  e,  e  real  and  positive  be  given.  There  exist  state 


functions  ©j)  e  S ,  defined  for  each  tuple 


maxi 

t 


j||©j('r))- Ja,|^,(T))|-dT 


<  £ 


(26) 


for  all  (a„  ...,  . 

The  proof  of  this  lemma,  while  not  difficult,  requires  extensive  manipulations.  We  will 
provide  it  in  a  companion  report  devoted  exclusively  to  the  chattering  aspects  of  our 
proposed  architecture.  Notice  that  the  lemma  says  that  every  chattering  combination  of 
state  functions  of  the  form  of  (17)  on  a  set  S,  can  be  realized  as  a  linear  combination  of 
elements  from  S  with  an  arbitrary  small  error  in  the  integral  sense  (see  (26)).  We  also 
have  the  fact  that  under  strong  continuity  assumptions  on  the  state  functions  the  converse 
of  the  lemma  is  also  true.  By  choosing  the  boundary  conditions  in  our  lattice  model 
appropriately,  this  assumption  is  not  limiting  . 

Now  we  proceed  to  describe  the  functionality  of  the  Input  selector  function.  This 
function  models  the  interaction  of  the  particle  with  the  excitation  and  with  the  other 
particles  in  the  lattice.  Some  examples  of  possible  peirticle  interaction  are  shown  in  Figure 
8.  For  simplicity,  only  nearest  neighboring  interactions  are  shown,  but  the  model  is  not 
limited  to  these  cases.  The  central  task  implemented  by  the  input  selector  function  is  to 
establish  on  one  hand  the  interaction  of  a  particle  with  the  other  particles  in  the  lattice  and 
on  the  other  hand  the  interaction  of  a  particle  with  the  excitation  field  in  order  to 
characterize  these  two  tasks  we  look  more  closely  at  the  excitation  field  Hamiltonian  Hj . 

Locally,  centered  in  a  particle,  it  is  convenient  to  write  Hj  as  the  sum  of  two  terms,  the 
first,  Hg,  corresponding  to  the  interaction  of  the  particle  with  the  excitation,  and  the 

second  H,  corresponding  to  the  interaction  of  the  particle  with  its  internal  energy  and 
with  the  interaction  energies  of  the  other  particles  in  the  lattice.  Thus  we  write 

Hf  =  Hg  +  H,  (27) 
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Figure  8.  Some  examples  of  particle  interaction 


Under  general  assumptions  about  the  lattice,  the  interaction  Hamiltonian  H,  for  given 
particle  can  be  written  as 


H,|>F)=  Hi-)|'P> 


(28) 


where  Hj^  for  each  j  is  the  interaction  Hamiltonian  of  the  given  particle  with  its 
'neighboring'  particle  ij  ,  and  g4,  is  the  neighborhood  function  at  the  current  state  of  the 
particle,  The  neighborhood  function  represents  the  local  (at  the  current  state 
structure  of  interaction  of  the  particles  in  the  lattice.  This  structure  can  be  determined  by 
building  the  lagrangian  associated  with  the  lattice  and  go  to  the  procedure  of  quantization 
that  we  discussed  earlier.  We  will  carry  out  this  task  once  the  details  of  the  physics  of  the 
lattice  are  defined  in  the  evaluation  phase.  The  approach  consist  in  defining  potentials  to 
characterize  interaction.  For  example,  a  dipole  (particle-to-particle)  interaction  between  a 
particle  located  in  coordinates  x,  y  of  the  lattice  may  be  characterized  by  a  potential  V  of 
the  form: 


where  d,  „  is  a  relaxation  function. 

A-y 

For  the  purposes  of  analysis  and  also  to  determine  a  detailed  formulation  of  the 
realization  step  in  our  paradigm,  it  is  convenient  to  assume  that  a  finite  set  of  primitive 

excitation  Hamiltonians  E  =  {H^  ,  i  =  1,...,  n^}  can  be  realized  and  that  implementation 
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of  our  excitation  is  carried  out  by  chattering  among  the  elements  of  E  over  the  update 
interval  A .  Specifically, 


H.I'Flf)) 


H;>|'P(T))Te 


(29) 


Where  the  sets  1;  (t)  are  defined  by  expression  (18).  This  type  of  probabilistic  resonance 

is  central  to  our  proposed  implementation  of  the  quantum  processor.  The  idea  is  to  induce 
a  probability  distribution  p  on  the  states  of  the  particles  on  the  lattice  so  that  the 
realization  criterion  is  satisfied. 

Thus  a  computation  in  the  lattice  is  a  propagating  probabilistic  wave-train  in  which  the 
state  of  each  particle  is  the  probability  distribution  of  its  pure  states.  This  is  illustrated  in 
Figure  9.  Specifically  the  lattice  is  at  an  initial  probabilistic  state,  the  programmed 
excitation  is  impinged  and  after  a  transient  period  ^  has  elapsed,  the  read  out  period,  the 
transducer  is  activated  to  effect  the  eigen-value  observation.  After  the  observation  has 
been  made,  the  computation  process  is  complete.  If  the  results  are  not  satisfactory  the 
computation  is  started  again  from  the  initial  probabilistic  state  and  the  read-out  period  is 

extended  to  ^,  >  ^  .  This  extension  period  cannot  be  extended  arbitrarily  because  the 
thermal  relaxation  mechanisms  in  the  lattice  will  induce  eventually  dejsoherence,  that  is, 
the  loss  of  the  probabilistic  resonance  described  above. 

Let  {Ej}  be  the  discrete  spectrum  of  the  physical  Hamiltonian  Hq.  Then  it  can  be  shown 
using  (24),  that  the  transition  probability  from  eigenstate  i  to  eigenstate)  satisfies, 

ih|pij(t)  =  (E,  -  Ej).  p,  +  (30) 

In  [1 1],  it  is  shown  that  if  the  system  is  at  state  at  time  t=  0,  then  the  presence  of  the 
relaxation  Hamiltonian  (28)  causes  the  corresponding  probability  density  term  p^^  to 
decay  exponentially  with  time.  For  small  values  of  t,  p-^  is  the  largest  term  in  the  matrix 

representation  of  the  probability  density  operator.  Assuming  no  excitation,  and  using  (30) 
this  term  satisfies  the  following  equation: 

-(>*(».)«)  01) 

with 

ih|p,  =  (Ej-  E,).  p,  -  Pj.(H,)„  fori* j  (32) 
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The  solution  of  (31)  and  (32)  for  p--  is  given  by: 

Pii(t)= 


(33) 


where  |  is  a  relaxation  time  that  depends  on  the  physical  characteristics  of  the  lattice. 
Thus  the  effect  of  the  inter  particle  excitation  Hamiltonian  is  to  randomize  the  state 
transitions  in  each  particle.  The  computation  described  above  will  be  corrupted  by  this 

effect.  Therefore,  the  read-out  period  ^  must  be  chosen  so  that  this  effect  is  acceptable; 

this  means  ^  We  will  use  this  bound  in  our  design  specifications  of  the  quantum 

processor  and  also,  in  our  computability  analysis. 


Forward  Wave  Direction 


Backwards  Wave  Direction 


Figure  9.  Wave  Propagation 
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Conclusions 


In  this  paper  we  present  a  preliminary  design  of  a  hardware  architecture  for  computing 
initial  segments  of  primitive  recursive  functions  and  iterative  processes.  The  formulation 
of  the  architecture  is  based  in  a  paradigm  which  proposes  a  procedure  for  1-  encoding  a 
function  or  process  and  2-  Carry  out  the  computation.  The  paradigm  is  firmly  rooted  in 
the  formalism  of  quantum  mechanics.  We  propose  as  our  representation  of  the 
architecture  a  generic  regular  multiparticle,  two-dimensional  lattice.  This  lattice  is  a 
model  of  crystal  structures  that  in  principle,  can  be  produced  in  the  lab  today. 
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